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LHCb strategies for 7 from 5 — > DK 

Y. Xie a (on behalf of the LHCb Collaboration) 

a School of Physics, University of Edinburgh, Mayfield Road, Edinburgh EH9 3JZ, UK 

One of the most promising ways to determine the angle 7 of the CKM unitarity triangle is through measurement 
of the tree level processes B — > DK. The LHCb collaboration has studied the potential of these decays employing 
the Atwood-Dunietz-Soni (ADS) and Dalitz methods, making use of a large sample of simulated data. For each 
method the expected sensitivities to the angle 7 are presented in this report. 
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1. INTRODUCTION 

This report is arranged as follows. In Section [1] 
we present a physics motivation, followed by a 
general discussion of B — > DK decays 0. Sec- 
tion[2]gives a brief description of the LHCb detec- 
tor, the simulation and the event selection tech- 
niques. An introduction to the ADS method to 
extract the CKM angle 7 from B — > DK decays 
and its application and expected performance in 
LHCb is presented in Section [3J Section 2] de- 
scribes the use of the Dalitz method to extract 7 
from B — > DK decays and the expected LHCb 
sensitivities. We conclude in Section [5l 

1.1. Motivation 

In B hadron decays, tree level processes are 
generally dominated by Standard Model contri- 
butions while new physics mainly affects loop di- 
agrams. Any difference between the CKM uni- 
tarity triangle as measured in tree level processes 
compared to loop level processes would indicate 
new physics in the flavour sector. Currently the 
unitarity triangle as determined by the tree level 
measurements of 7 and the CKM matrix element 
\V u b\ is consistent with the triangle determined 
by measurements of the loop dominated parame- 
ters efe, Arrid, Am s and sin 2/3 pQ. However, the 
present measurements of 7 from B — > DK at B 
factories have large errors: 7 = (92 ± 41 ± 11 ± 
12)°(BaBar) 0,7= (53tJ| ± 3 ± 9)° (Belle) [3J. 

The precision on the angle 7 from B — > DK has 

1 In this report D represents a D° or D°, B represents a 
B ± , B° or B° and K represents a K ± , K*° or K*°. 



to be improved to at least 5° p] to match the pre- 
cision of its indirect estimate from a global fit to 
CKM parameters excluding direct measurements 
of 7. 

In LHCb, the tree level decays B° s -> DfK T 
can also be used to determine the angle 7 in a 
theoretically clean way. This measurement is ex- 
pected to have a sensitivity of 14° with 2 fb _1 of 
data. We will see below that the B — > DK decays 
have a greater statistical sensitivity to the angle 
7 with an equivalent amount of data. 



1.2. Features of B — > DK decays 

The CKM favoured process B -> D°K and 
disfavoured process B —> D°K can be described 
with three parameters: a weak phase difference 
7 = arg{-{V ud V* b / {V cd V* b )), a strong phase dif- 
ference 5b and the ratio of magnitudes between 
the disfavoured and favoured amplitudes, defined 
as rs- If the D° and D° decay to a common fi- 
nal state, then the interference between the two 
amplitudes via D° and D° allows the extraction 
the angle 7 with several methods, as illustrated 
below. 

It is expected that rs is small in the B ± — > 
DK^ case and has a value of about 0.1 [415] due 
to colour suppression in the CKM disfavoured 
amplitude. For B° — » DK*°, both amplitudes 
are colour suppressed, therefore tb is expected to 
be larger [S]. 
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2. THE LHCb DETECTOR 

The LHCb detector is a single arm spectrom- 
eter dedicated to the study of CP violation in 
B meson decays at the Large Hadron Collider, 
which will start operation at CERN in 2007. 
The detector and its expected performance is de- 
scribed in detail in [6]. Here we only empha- 
size that LHCb has a 94% tracking efficiency for 
tracks with momentum above 10 GeV/c, a 93% 
identification efficiency and a corresponding 
probabilty of 4.7% for a to be misidentified as 
a K ± for the momentum range 2 — 100 GeV/c. 

2.1. Data simulation and event selection 

Monte Carlo simulation data produced with 
Pythia and Geant4 are used to study the trig- 
ger, the reconstruction and event selection, which 
in turn allows the physics performance to be as- 
sessed. 

We use event samples consisting of 260 mil- 
lion minimum bias event for trigger studies, 140 
million inclusive bb events and dedicated signal 
events for selection and background studies. 

Sensitivities on physics parameters are ob- 
tained using fast simulation. These are based on 
efficiencies, resolutions and background levels ob- 
tained from the full simulation. 

In our studies we use the following information 
to discriminate between signal and background 
events: charged particle identification informa- 
tion based on the Ring Imaging Cherenkov detec- 
tors; invariant masses; impact parameters; trans- 
verse momenta; x 2 of decay vertices of the B, D, 
K* and Kg particles; the opening angle between 
the momentum direction of a B and its flight di- 
rection; the event topology itself. 

All selection cuts are optimized to reject most 
background events from a large sample of inclu- 
sive bb events while retaining a reasonable sig- 
nal efficiency. The ratio of background to signal, 
B/S, is assessed using a separate inclusive bb sam- 
ple. 

3. ADS METHOD AND SENSITIVITY 

3.1. Description of the method 

Atwood, Dunietz and Soni (ADS) [7 suggested 
a method of determining 7 based on the recon- 



struction of non-CP eigenstates for decays com- 
mon to both D° and D°. An example is the 
hadronic final state Kit, which may arise from 
the Cabibbo favoured (CF) decay D° -> K + Tr~ 
or the doubly-Cabibbo suppressed (DCS) decay 
D° — > K + ir~ . The relation between the CF and 
the DCS decay is described by a magnitude ratio 
r^ 7r and a strong phase difference 5^ n : 

K i5 K, A(D° -» K-n+) = A(D° -> K+n-) 
D A(D°^K-tt+) A(D° -> K+n-)' [ ' 

Therefore there are four possible B decays, 
whose decay rates can be written as follows 

T(B~ -» (K-tt+) d K-) cx 
1 + {r B r%*f + 1r B r$ w cos(<5 B - 8%" - 7), (2) 

T(B~ -> {K + tt') d K~) cx 
r% + {r^f + 2r B r^ cos(<5 B + 6%" - 7), (3) 

T(B+ -> {K + n-) D K + ) cx 

1 + (r B r^) 2 + 2r fl r£" cos(<5 B - 5%" + 7), (4) 

T(B+ -» (K-tt + ) d K+) cx 

r| + (r£*) a + 2r B r« * cos(6 B + 5%* + 7), (5) 

where the constant of proportionality is the same 
in each expression. The relative rates of the four 
processes yield three observables which depend 
on five parameters 7, r B , r|j , 5 B and , of 
which rjj is already well known. It is necessary 
to use different D decays in order to determine 
all parameters. 

Similarly, the four-body decay D — » Ktttttt 
provides three new observables which depend on 
7, r B , 8b and two new parameters r^ 37r and 
§K3tt |2| Further information can be added by 
including D decays to CP-eigenstates, such as 
K + K~ and Tt + n~. Each provides one more ob- 
servable without introducing any new parame- 
ters: 

T(B~ -> (h+h-) D K-) cx \+{r B f+2r B cos(5 B -7), (6) 
T{B+ -> {h+h-) D K+) cx l+(r B ) 2 +2r B cos(5 B +j). (7) 

2 In fact the inclusion of D° — > Kmrn will bring many 
more parameters - and additional information - because of 
the resonant substructure of the decay. This complication 
is neglected in the present discussion. 
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3.2. Performance with charged B mesons 

Results from the B factories favour a small 
value of r B for charged B decays [819]. We set 
tb = 0.077 for this study. The following assump- 
tions are made: = r£ 3n = 0.06, -25° < 
< 25° and -180° < 6% 37T < 180°. 

Tab. Q] shows the expected signal and back- 
ground yields in 2 fb _1 of data. Of the 17.7 k 
background events in the favoured B — > (Ktt) dK 
modes, 17.0 k are from the decay B — > Dir, 
which has a 13 times larger branching ratio, with 
a 7r misidentified as a K, and 0.7 k are combi- 
natorial background events. In the suppressed 
B — » (Ktt)jjK modes, the combinatorial back- 
ground dominates. 

Our simulation study indicates that the signal 
yields and background level are very similar for 
the D — > K3ir modes. Therefore we use the same 
yields and B/S as in D — > i^7r. 

A large number of fast samples are generated 
based on the yields and background levels given 
above to estimate the statistical precision of 7. 
As shown in Tab. [2] a precision of 5° — 15° for 
7 is achievable for 2 fb _1 of data, depending on 
the parameter values of 8^ and S^ 37T and with 
rs = 0.077. Better precision can be achieved for 
larger re values. 

LHCb is also investigating the feasibility to in- 
clude the decay B ± -> D*°K ± in the ADS analy- 
sis. An attractive feature of the D*° meson is that 
it can decay to two final states D°7r° and D°j. 
These have opposite CP eigenvalues, which lead 
to a difference of it in their strong phases. This is 
a useful constraint if the two decays can be dis- 
tinguished experimentally. However, these decays 
are difficult to fully reconstruct in LHCb because 
the detection efficiency of a soft photon in the 
electromagnetic calorimeter is very low while the 
background is enormous. Alternative approaches 
which employ a partial reconstruction or which 
make use of the event topology to reconstruct the 
momentum of the ir° /j are under study. 

3.3. Performance with neutral B mesons 

The same method can also be applied to the 
decay B° -> DK*° with D -> Ktt, KK or tttt. 
Assuming r B = 0.4, 55° < 7 < 105° and -20° < 
5b < 20°, the expected signal yields in 2 fb _1 



of data and the background-to-signal ratio B/S 
are given in Tab.[3l The corresponding statistical 
precision of the angle 7 is expected to be 7° — 10°. 

4. DALITZ METHOD AND SENSITIV- 
ITY 

4.1. Three-body D decay 

This method for determining 7 was proposed 
in [TU]. It makes use of the decay B ± — ► DK ± 
followed by a multibody D decay into a CP 
eigenstate Here we explain the basic idea us- 
ing D — > Kgir + ir~ as an example. The Dalitz 
phase space of this decay D — > ATiy7r + 7r~ can 
be fully parameterized with two effective masses 
m\ = m 2 (K° s TT+) and m 2 _ = m 2 (K%ir-). The 
D° and D° decay amplitudes can be written as 
functions f{m_,rnA_) and f(m 2 lr ,m?_). 

The total Dalitz decay amplitudes, defined as 
A- = A{B- -> (K° s -k+-k ) D K~) and A+ = 
A(B + (Kgir + ir~)B.K + ), are sums of contri- 
butions via D° and D°: 

A- = /(ml^rnl) +r B e l ^ +5 ^f(ml,m 2 _), (8) 
A+ = f(m 2 + ,m 2 _) + r B e l{l+5B) f(m 2 _,m 2 + ). (9) 

In the isobar model [TT] f(m 2 +: m 2 _) is a coher- 
ent sum of contributions of different resonances: 

N 

f(m 2 + ,m 2 _) = ^2a je za i Ajimltml) + be tfJ 7 (10) 
3=1 

where a,j , otj, b and (3 are model parameters which 
have been measured well at the B factories [213] . 
The B^ decay rates are given by 

T~(m 2 + ,m 2 _) = \f(m 2 _, m 2 + )\ 2 +r 2 B \.f(m 2 + ,m 2 _)\ 2 

+ 2rBRe[f\m 2 _,m 2 + )f(m 2 + ,m 2 _)e i{ ~ 1+&B \ (11) 

r + (m 2 + ,m 2 _) = \f(m 2 + ,m 2 _)\ 2 +r 2 B \f(m 2 _,m 2 + )\ 2 

+ 2r B Re[f*(m 2 + ,m 2 _)f(m 2 _,m 2 + )e i ^ +SB) }. (12) 

We can see that the interference terms are sen- 
sitive to 7, which can be determined by measuring 
r~ (m 2 + ,m 2 _) and T + {m 2 + ,m 2 _) across the Dalitz 
phase space. 
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Table 1 

Expected signal yields S, number of background events B and the rato B/S in 2 fb _1 of data for ADS 
decay modes of B^ corresponding to 8 b = 130° and 8^ — 0°. The uncertainty on the background 
estimates is around 60% for the rarest modes. 



decay mode 


5 


B 


B/S 


B+ -> (K+tt-) d K+ 


28 fc 


17.7 fc 


0.6 


B~ -» (K-tt+) d K- 


28 fc 


17.7 fc 


0.6 


B+ -> {K-it+) D K+ 


530 


770 


1.5 


B~ -> (K+n-) D K- 


180 


770 


4.3 


B+ -> (X+is:-/7r+7r-)Dif + 


4.3 fc 


4.3 k 


1.0 




3.3 fc 


3.3 k 


1.0 



Table 2 

The statistical error of 7 for different values of 8^ and S^ 371 for = 0.077. Numbers with * are RMS 
values quoted for non-Gaussian distribution of fit results due to close lying ambiguous solutions. These 



will disappear 


as the si{ 


mal yields 


increase. 










sKtt 

°d 


-25° 


-16.6° 


-8.3° 


0° 


8.3° 


16.6° 


25° 


8^ = -180° 


8.6° 


7.5° 


6.5° 


6.8°* 


7.2°* 


7.3°* 


6.0°* 


8^ = -120° 


6.0° 


6.3° 


6.3° 


6.4° 


6.2° 


6.2° 


4.7° 


5 K3tt = _g o 


8.0° 


7.9° 


8.1° 


7.8° 


7.4° 


6.7° 


6.2° 


5% 371 = 0° 


10.3°* 


11.1°* 


12.°* 


11.5°* 


12.1°* 


13.1°* 


13.0°* 


8% 37T = 60° 


9.1° 


10.6° 


11.2° 


12.9° 


13.4°* 


15.0°* 


15.2°* 


5% 37r = 120° 


11.6°* 


11.3°* 


11.8°* 


11.0°* 


10.9°* 


11.1° 


10.9° 


5§ 3 * = 180° 


8.5° 


7.4° 


6.5° 


6.8° 


7.1° 


7.3° 


6.5° 



To reconstruct the B ± -» (K%-k + -k~)dK ± 
events is challenging with the LHCb detector as 
only 25% of the Kg particles decay inside the 
active region of the vertex detector. The ex- 
pected signal yield in 2 fb _1 of data will vary 
between 1.5 k and 5 k, depending on how many 
of the K s decays successfully found offline can 
be reconstructed within the CPU constraints of 
the High Level Trigger. A full simulation has 
been performed to estimate the background level. 
The combinatorial background is expected to con- 
tribute less than 3.5 k events and contamina- 
tion from B^ — > (i^^7r + 7r _ )i)7r ± is expected to 
be around 1200 events. Our present sensitivity 
studies for 7 do not include background, and do 
not take into account the non-flat acceptance ef- 
ficiency in the Dalitz space, which is expected as 
a result of the trigger and offline selection. Un- 
der these assumptions, a statistical precision of 



er 7 w 8° - 16° is achievable in 2 ffcT 1 of data. The 
actual statistical precision will depend on the fi- 
nal background level and on re- 

LHCb is also investigating the decays B ± — > 
{K S K+K-) D K ± and B° -> (K s ir+w-) D K *°. 

4.2. Four-body D decay 

The Dalitz method can be extended from three 
to four-body D decays. In this case five param- 
eters are required to describe the Dalitz phase 
space. The D decay model has been studied in 
the FOCUS experiment [12] and our 7 sensitivity 
studies are based on their results. 

Assuming a branching ratio of B(B ± — > 
(K + K-tt+tt-) d K ± ) = 9.5 x 10~ 7 , our full simu- 
lation yields 1.7 k events in 2 fb" 1 of data. Based 
on this we estimate a statistical precision for 7 to 
be <r 7 « 14°, where we have assumed 7 = 60°, 
tb = 0.08 and 8b = 130°, and have not yet in- 
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Table 3 

Expected signal yields in 2 fb _1 of data and the background-to-signal ratio B/S for the ADS modes of 
neutral B mesons. Upper limits with 90% confidence level are quoted for B/S. 



decay mode 




S 


B/S 


B° -> {K-n+ 


) D K*° + c.c. 


3400 


<0.3 


B° -> {K+tt- 


) D K*° + c.c. 


500 


<1.7 


B° -> {K+K- 


' /ir + ir-) D K* a + c.c. 


500 


<1.4 



eluded background and detector effects. Results 
from more recent studies can be found in [13] . 

We are also studying an amplitude analysis for 
B ± -> (K ± tt t tt + tt~) d K ± . Compared with the 
ADS analysis discussed in Section [3j the advan- 
tage of this method is that it takes into account 
the variation of the strong phase <5^' 37r in the 
Dalitz phase space. 

4.3. Systematic errors 

The biggest systematic uncertainty in the 
Dalitz method arise from the model depen- 
dence of the D decay. In the present B-factory 
£? ± — » (Kgir + ir~)£)K ± analyses this uncertainty 
is around 10°. This error is expected to reduce 
significantly through exploitation of the coher- 
ently produced D mesons available at CLEO- 
c [14] and BES [15] . A discussion of how these 
data may be used in a model independent analysis 
can be found in [10116] . Similar techniques can be 
used for the four-body decay mode, where LHCb 
also expects large numbers of flavour-tagged D 
decays for use in model calibration. 

5. CONCLUSIONS 

We have shown that LHCb will be able to 
extract the CKM angle 7 in several ways with 
B — > DK decays. The combined result is ex- 
pected to have a precision of around 5° with 
2 fb _1 of data. Such a result will make it pos- 
sible to compare the LHCb measurement of the 
angle 7 with the indirect determination from a 
CKM fit and thereby perform a stringent test of 
the Standard Model. Together with improvement 
on the \V u b\ measurement at the B factories this 
will provide a precise reference Unitarity Triangle 
against which new physics searches can be com- 
pared. 
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